Europhysics Letters 



PREPRINT 



Droplet rebound from rigid or strongly hydrophobic sur- 
faces 

A. GOPINATH 

Division of Engineering and Applied Sciences, Harvard University, Cambridge, MA 
02138 



PACS. nn.mm.xx - First pacs description. 
PACS. nn.mm.xx - Second pacs description. 
PACS. nn.mm.xx - Third pacs description. 



Abstract. - A weakly deformable droplet impinging on a rigid surface rebounds if the surface 
is intrinsically hydrophobic or if the gas film trapped underneath the droplet is able to keep the 
interfaces from touching. A simple, physically motivated model inspired by analysis of droplets 
colliding with deformable interfaces is proposed in order to investigate the dynamics of the 
rebound process and the effects of gravity. The analysis yields estimates of the bounce time 
that are in very good quantitative agreement with recent experimental data (Okumura et. al., 
(2003)) and provides significant improvement over simple scaling results. 



Introduction. - The collision and rebound of droplets has been a subject of interest and 
investigation for many decades. Impacting drops undergo energetic bounces if the pressure 
in the gas film separating the drops deforms the drop surfaces sufficiently to transform the 
drop's kinetic energy into deformation energy before contact occurs. Collisions in a incom- 
pressible, continuum gas usually result in a bounce with contact occurring due to inter-particle 
forces or surface imperfections^'^. Bouncing can also result when drops impinge on rigid super 
hydrophobic'^ surfaces or on rigid walls heated to well above the Leidenfrost temperature*'^. 
Very recently, Okumura et. al.'^, presented experimental data for the collision and rebound, 
in air, of 400 — 1000 /im water drops from a super-hydrophobic surface. They found that the 
measured bounce time was consistently larger than that predicted by simple scaling. Further- 
more, with the drop radius held fixed, the contact time was seen to increase significantly as 
the velocity of impact decreased. This was attributed to effects of finite drop size and gravity. 

In this letter, a physically motivated model is proposed to investigate this scenario, specif- 
ically the dependence of the bounce time and deformation on droplet size and velocity. Use of 
the methodology developed in similar problems involving droplets colliding with deformable 
interfaces"'^'^, indicates that even when gravity is absent the bounce time for a liquid drop 
increases as the impact velocity decreases. Consideration of gravity induced effects suggests 
that the bounce time is modestly modified relative to the zero-gravity value. The predicted 
results are in excellent quantitative agreement with experimental data of Okumura et. al.,^ 
showing that both effects are needed to account properly for the details of the bounce. 
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Fig. 1 - (a) Schematic of a drop of initial radius a approaching a rigid hydrophobic surface with scaled 
velocity Wd- (b) A close up of the inner region reveals a dimpled drop with a clearly defined rim. For 
a drop supported by a film of gas, the gas pressure is 0(2a/a). 



Dimensionless analysis yields reduced parameter space. - Consider a droplet of radius 
a comprised of liquid with density pd, viscosity pd and interfacial tension cr colliding with a 
rigid flat surface at impact speed Uc- Let the ambient gas be incompressible and continuum 
with density pg, pressure Po and viscosity pg. Dimensional analysis then indicates that the 
parameters determining the characteristics of the rebound process are - (a) the Weber numbers 
based on drop and gas properties, Wd = pdU^acr and Wg = pgll^a/a that measure the relative 
importance of inertial pressures and surface pressure, (b) the capillary numbers based on 
drop and gas viscosities, Cd = PdUc/'^ and Cd = PgUc/(J, and (c) the Reynolds numbers 
Rg = aUcPg/pg and Rd = aUcPd/pd characterizing viscous effects in the gas and inside the 
drop. The effects of gravity are quantified by the Froude number, Fd = U^/{ag). Droplet 
inertia is quantified by the Stokes number S ^ 0{Wd/Cg),tha,t also measures the effectiveness 
of gas viscosity in dissipating the kinetic energy of the drop. 

A large number of physically relevant collisions occur at conditions satisfying Wd <IC 1, 

1/2 

5 > 1, i?d ^ 0, i?g < 1, W^' /Rd and Cg < 1. In such cases, viscous dissipation inside 
the drop and in the gas outside may be neglected during the collision process as long as the 
droplet rebounds^. That is, the dissipation in the gas film has a small impact on the bounce 
provided we are far from the bounce-contact transition. The primary role of the film is then 
to cushion the droplet from physical contact with the surface. We expect the dynamics to 
primarily depend on just two parameters - Wd and Fd and seek a theoretical description of 
the approach and rebound in this special limit. 

Collisions in zero gravity. - Consider for now collisions in gravity-less conditions so that 
Fd = oo. It is useful to first visualize the approach and rebound scenario - for simplicity, let us 
choose a collision occurring in the presence of a gas. For 5^1, changes in the droplet velocity 
occur only when the gap thickness, ho <C a, as a result of the large lubrication force in the gap. 

1/2 

Drop deformation becomes important when the gap thickness is 0{aCg ) <^ a. This is easily 
seen by balancing the viscous force exerted by the gas with the capillary pressure 0{2a/a). 
The center of mass speed is chosen to be Uc at this point. As the gap becomes smaller, the 
confined air film leads to a strong localized pressure that is 0{2a / a), that when seen from far 
seems to act over a temporally varying nearly flat region of radial extent rd{t) as depicted in 
Figure (la). Closer examination of the near-contact region as depicted in Figure (lb) reveals 
that the inner region may in fact have a more complicated shape - for instance during the 
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approach, a dimple like shape is seen with a radius of roughly rd{t) and a gap thickness that 
is a maximum at the centerline. Insofar as the bounce process itself is concerned, the detailed 
shape of the inner region is of secondary importance - what is important is the angular extent 
of the flat region^ As the radius of the flat region increases, the center of mass moves slower 
as as more of the initial kinetic energy is channeled into surface deformation modes. At some 
point, the centers of mass of the drops reverse their direction of motion and the drops begin 
to rebound - during this process the flat region begins to shrink and finally vanishes. For a 
droplet colliding with a super-hydrophobic surface we expect a similar picture to apply. The 
reason being that due to the nature of the surface-drop interaction, the contact angle is nearly 
180° and the near-contact area is flat to leading order here as well. Thus, drop is nearly 
spherical except at a small region where it is flat. This shape is reminiscent of the shape 
found by Mahadevan and Pomeau* for small drops rolling due to gravity. 

Simple scalings^'^'^ for the bounce time and deformation are readily obtained. For the 
deformation energy to balance the 0{2Trpda^U^ /3) initial kinetic energy, the characteristic 
magnitude of the interface deformation should scale as 0{aWy^) ^ 0{pdU'^a'^ / crY^'^ ■ The 
time scale for the deformation to relax is approximately 0{aWy^U~^) ~ 0{pdO? j 0)^!"^ . 
This value, which we denote by tr, was flrst obtained by Lord Rayleigh in his study of the 
vibrational modes of a drop^. Note that the Rayleigh time scale is independent of the impact 
velocity. The angular extent of the flat region is obtained by a force balance on the drop and 
is 0(aWel^). 

To obtain a more accurate picture of the rebound process, the equations for the flow inside 
the drop have to be solved. The formulation follows directly from earlier work on collisions 
with deformable interfaces''^. A spherical coordinate system that is fixed in space and time 
in the interior of the droplet, such that the origin lies at a distance a from the rigid surface, 
is chosen. The velocity field in the drop is divergence free and governed by the Navier-Stokes 
equation. Scaling all lengths by a, time with aW e^j'^U'^^ , and the pressure by {pdU^W^ ^^^) 
and recognizing that the 0{a? pdP^^) time scale for viscous diffusion of momentum in the drop 
is much larger than the bounce time, we find that the fluid in the drop undergoes inviscid 
impulsive motion described by, 

9tUd = -Vpd. (1) 

Invoking the concept of a velocity potential, <f) such that Ud = V(^, we find that = and 
the dimensional pressure field is given by 

p'd=p'o + '^{o/a)-pddtcj>' 

being the ambient pressure. The effect of the 0{2a/a) pressure field localized in a flat 
region of angular extent a{t) is incorporated via the equation p'g = p'^-\- 2(a/a)H(a(f) — Q), 
9 being the angular co-ordinate of a cylindrical co-ordinate system as in Figure (la), and H 
being the Heaviside function. The radial position vector of a surface point at r; = cos 6 and 
the velocity potential are expanded using the Legendre functions Pk of order k and are given 

by 

00 00 

r:{r],t)=aR{r],t) = a{l + W^/^Y.^k{t)Pk{v)) and <^ = ^ Bfc(i)(r*)'=Pfc. (2) 

fe=0 fc=0 

Application of the kinematic and normal stress boundary conditions, the constant volume 
constraint along with simple geometry then yields 

dt{Bk) = -[Pfc_i(cosa) - Vk+i{cos a)]W-'^^ - [k{k + 1) - 2]Dk, dt^D^) = kBk, 
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Fig. 2 - Scaled rebound time, Tj = TI/th as a function of Weber number, Wd- The solid curve is 
the theoretical prediction for S — > oo, Wd ^ 1 and Fd = oo. The symbols are experimental data 
points'' for water drops in air. The squares correspond to a = 0.6 mm, the stars to o = 0.4 mm and 
the circles to a = 1 mm. 



and 1 + ^ Dk{t)Vk{cosa) = (cosa)-^ (3) 
fc=i 

Equation set (3) is then solved in conjunction with the initial conditions Dk>i{t = 0) = 0, 
Bk>2{t = 0) = and Bi{t — Q) — 1. The rebound is complete when the center of mass of 
the droplet reaches its initial position viz., Di{t = Tb) = 0. Numerical solutions are obtained 
for various values of Wd by using a second-order predictor-corrector method and varying the 
number of modes until convergent results are obtained. 

Figure (2) summarizes the results for — oo and compares them to experimental obser- 
vations of Okumura et.al.^. The solid line corresponds to the theoretical prediction obtained 
by numerically solving equation (3). We notice immediately that in spite of the different 
impact velocities and sizes, there is a reasonable collapse of the experimental data around the 
theoretical prediction, indicating that Wd is the relevant parameter to use. Good qualitative 
and quantitative agreement is seen, especially for the smaller size drops. Some of the data 
points for the 1 mm drop do differ considerably from the predicted values but the reasons for 
these are not clear. Nonetheless we find that the theoretical result yields a good prediction 
even as the Weber number becomes 0(0. 1). 

Both experimental data (that include gravity effects) and the numerical results (with- 
out gravity) consistently indicate that as Wd decreases, the bounce time increases with 
decreasing Uc- To investigate this behavior in more detail, analytical asymptotes for the 
bounce time were obtained by a singular perturbation analysis of (3). It is found that when 
la}^"^ (W^ ^^'^) ^ 1, all the deformation modes do not scale similarly as 0(aW^''^), instead 
a separation of scales develops. The fc = 1 mode corresponding to translation of the center 
of mass scales as 0(aW^^^ In^^^ (W^ ^^^)) while the fc > 2 surface deformation modes scale 
as 0{aWy^ hi~^^^ {W^ ^^^))- As a result, when Wd — * 0, the highly localized pressure field 
causes the k > 2 flow modes to become negligibly small. The radial extent of the inner region 
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then couples directly to the A; = 1 translational mode. In other words, rather than deform 
uniformly, the approaching drop prefers to flatten the near contact region and shift the center 
of mass appropriately in response to the ensuing force. A spring like motion manifests itself 
wherein the center of mass executes a harmonic oscillation^'^ but like a non-linear spring that 
stiffens strongly as Wd O.In fact, for asymptotically small Wd, the bounce time scales as, 

n « 7^^/(273) [ aU-'wy' In^^ (i^" V4) j (4) 

while the maximum extent of the flat region scales as amax = [ 2Wd/(3 In ]^/^. As 

the Weber number increases, all the deformation modes become equally important and the 
separation of scales diminshcs. The asymptotic expression in (4) consistently undcr-prcdicts 
the actual bounce time since the analytical approximation is obtained by ignoring the k = 2 
modes. For finite Wd, the symmetry between approach and rebound is broken. The effects 
of finite Weber number render the collision process asymmetric, increasingly so as the Weber 
number approaches 0(1) values. 

First effects of gravity. - The effects of gravity cannot be ignored when Wd ~ 0{F^) i.e, 
when the impact velocity is 0{c? pg^ j uY^"^ or smaller. To anticipate the effects of gravity we 
consider the case of weak but non-zero gravity. That is, Fd is finite but e\ = W^ F^ -C 
1. Thus the droplet is still weakly deformed and the flow and deformation modes are still 
described by (2). Note that the small parameter e\ is the ratio of the critical speed at which 
gravity becomes important and the impact speed - that is e\ = ag{apd/ cf)^^'^U~^ = Vg/Uc- A 
measure of the importance of gravity in affecting the bounce time is then obtained by scaling 
the actual bounce time, T^' with the gravity-free value and then plotting this normalized 
rcboimd time as a function of a suitable small parameter. Since the data for 400 micron 
drops correspond most closely to the limits we are investigating, these are chosen to test the 
theoretical predictions. 

Okumura et. al., proposed an approximate expression to incorporate the effects of gravity 
which may be written as 

(|,). = 2(l-icos-(^j^)). (5) 

In their expression, the zero-gravity bounce time, (T^°)o = {^.Str) is independent of Wd and 
thus does not depend on the impact velocity. The bounce time as predicted increases by at 
most a factor of 2 compared to the zero-gravity value. This expression is plotted as a dashed 
curve in figure (3a) and compared to the experimental obtained data points normalized the 
same way. The circles indeed map onto to small si values but the normalized bounce times 
are all much greater than unity. At first sight this may seem to imply that gravity plays a 
crucial role - however this is not quite correct. Since the true zero-gravity bounce time is 
dependent on Wd, a more accurate reflection of gravity effects is obtained upon incorporating 
this dependence. 

Figure 3b illustrates precisely this effect. The squares are the re-normalized data points 
obtained by calculating corresponding to the solid line in figure 2, and plotting T^/T^ as 

a function of the small parameter £2 = ^d^"^ In^^"^ i^d ^^^)-^d^ ■ d'^ta points map to low 
values of £2 as expected indicating the weak effects of gravity. More importantly, the scaled 
bounce times are within 20 % of the gravity-free value suggesting that the quantitative effects 
of gravity arc relatively modest. 

An approximate theoretical treatment allows us to anticipate the dependence of the scaled 
bounce time on £2 observed in figure 3b. For asymptotically small Weber number, the center 
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Fig. 3 - The effects of gravity on the rebound time are explored in this plot, (a) Comparison of 
equation (5) proposed by Okumura et. al. (dashed curve) with the normalized bounce times (open 
circles) for 400 micron droplets colliding in air. The x axis is £i = VgU^^ ■ The gravity-free bounce 
time is independent of Wd- (b) The solid line corresponds to equation (8) and incorporates the effect 
of Wd on the gravity-free bounce time. The squares are the same data points re-normalized differently 
and plotted as a function of £2 = Vg hi^^^ {W^^^'^)U^^ . 

of mass of the drop behaves as if attached to a spring with a Wd dependent stiffness. When 
the effects of an external gravity field are included, an additional force acts on the center 
of mass due to the weight of the drop. Set r*{ri,t) ~ a(l + SDiPkit]) + 0(62)) where = 
W]/^ ln^/2 (W^d"^^"*) and 62 = W^^ In'^/^ (W^cT^^"^)- The force balance on the center of mass 
of the drop has the form 

{ATTpday3){aSi)[U~\a5i) r^d^Di) = ~2{<j/a)TT{a6l/^)^a^ + {4TTpda^g/3){g ■ U,) (gt/,)"!. 

Algebraic manipulations using asymptotic properties of the Legendre polynomials show that 
a^{t) w Di{t) + 0(62) ■ The deformation mode Di then satisfies 

d',{D,) ^-^D, + -^(g . U,)(.gC/e)-^ = -^D, + 82 (6) 
ltd I 

the small parameter £2 = 5i/Fd and gravity has been chosen to act in the direction of the 
initial droplet motion. Equation (5) indicates that the relevant small parameter quantifying 
effects of gravity when Wd ^ 1 is not ei but €2. Solving (5) with the initial conditions 
Di{t = Q) = and {dtDi){t = 0) = 1 yields the solution 

Di{t) = l^e2+ [^(l + ^e^)]^sin(y|t + tan-i(-£2y|)). (7) 

This corresponds to the drop being spherical and moving with speed Uc when the collision 
process starts. 

It is easy to show from equation (6) that when < 62 ^ 1, the modified bounce time is 
given by 

^ = (l-|tan-i(-e2y|)). (8) 

It is clear that the bounce time increases relative to the zero-gravity value by a factor that is 
between 1 and 2, the latter value corresponding to £2 ^ 00. 
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Equation (8) shows surprisingly good agreement with the experimental data points which 
seems to suggest that gravity acts as a weak perturbation for these droplets. Perhaps, a more 
careful calculation in lieu of the simplified treatment just described will result in a better 
prediction for larger values of £2. Also note that the dashed and solid curves seem to have the 
same form - this is of course primarily due to the way we have chosen to normalize the curves. 
It should also be noted that if the direction in which gravity acts is not the same as that of 
the initial droplet motion, axi-symmetry is broken and the bounce times can differ noticeably. 

Summary and conclusions. - To summarize, the collision and rebound dynamics of 
weakly deforming drops impacting with velocities, Uc, satisfying (a^pg'^/a) <€. <€. {ap/cr)~^, 
is controlled by the interplay between the kinetic energy of the drop and energy stored in sur- 
face deformation due to finite surface tension, a with gravity acting as a weak perturbation. In 
the limit of asymptotically small Wd, the scaled bounce time T;, = Tl^{pa^ /a)~^/'^, diverges as 
In^^^ [ {pU'^a/ a)^^/'^ ] when Uc 0. To leading order, gravity tends to increase the time scale 
for the bounce relative to the zero-gravity value. Results of numerical calculations indicate 
that the bounce time increases relative to the asymptote as {ap/a)~^ - a manifesta- 

tion of the increasing effect of non-linear deformation modes and of the drop being deformed 
everywhere rather than merely locally. 

Is the 0(aWj^^) local flattening seen in the Wd "C 1 case, a precursor of the global 
deformation seen in pancake shaped drops for Wd > 1? It is tempting to think that as Wd 1, 
the maximum extent of the inner region grows till it becomes comparable to the drop radius a 
with the transition between the inner and outer region becoming gentler. Consider increasing 
the impact speed, Uc such that it becomes 0{ap/a)~^/'^ as happens when Wd 0(1). In this 
limit, gravity plays a negligible role. The strongly deforming droplet changes to a pancake 
like shape with characteristic radius Rd> a being different from the height Hd ^ 0{a^R^^). 
The initial kinetic energy is transferred mainly into internal flow modes within the drop. The 
maximum extent of the pancake drop is determined by a balance between the flow driven 
by inertial acceleration and the flow due to the surface tension induced pressure gradient 
at the edge of the pancake. This yields the balance {pdU'^/a) ^ {aR\/a^) thus yielding 
{Rd/ay ~ 0{pdU^a/a) ~ 0{Wd)- Although the scalings for the maximal apparent contact 
area seem to arise differently for cases Wd "C 1 and Wd ^ 1, it is interesting and in some 
sense satisfying to note that they exhibit the same dependence on Wd- 

* * * 

I would like to thank Prof. L. Mahadevan for his encouragement and help during the 
preparation of this manuscript. 

REFERENCES 

[1] Gopinath A and Koch D. L., J. Fluid. Mech., 454 (2002) 145. 
[2] Gopinath A and Koch D. L., Phys. Fluids, 13 (2001) 3526. 

[3] Okumura K., Chevy F., Richard D., Quere D. and Clanet C, Europhys. Lett., 62 (2003) 
237. 

[4] Karl A. and Frohn A., Phys. Fluids, 12 (2000) 785. 

[5] Yao S.-C. and Cai K. Y., Exp. Therm. Fluid Sci., 1 (1988) 363. 

[6] Rayleigh, Lord, Proc. R. Soc. Land. A., 29 (1879) 71. 

[7] Mahadevan L. and Pomeau Y., Phys. Fluids, 11 (1999) 2449. 



